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ABSTRACT 


A model is presented in which one calculates the paramagnetic resonance 
(epr) spectrum of vanadyl acetylacetonate (VAAC) dissolved in either a liquid 
crystal or isotropic solvent . It employs density matrix formulation in the 
rotating reference frame* The molecules occupy several discrete angles with., 
respect to the magnetic field and can relax to neighboring positions in a 
characteristic time x(0)o The form of t( 9) is found from a diffusion 
approach, and the magnitude of x(0) is a measure of how freely the VAAC probe 
tumbles in the solvent. Spectra are predicted for values of t between 
-11 -7 

10 sec and 10 sec. The epr spectrum, in the isotropic case, is obtained 
by summing the contributions from the allowed angles weighted by the polar 
volume element, sin 0. When applying the model to the nematic liquid crystal 
case it is also necessary to multiply by the Saupe distribution function. For 
this case, x(0) is obtained from the diffusion approach in which two diffusion 
constants are employed to reflect the difference in the parallel and perpendi- 
cular components of the viscosity. An attempt will be made to correlate the 
diffusion constant with the viscosity via the Einstein relation. 


Electron paramagnetic resonance (epr) has been used to measure the 
ordering of the probe molecule vanadyl acetylacetonate (VAAC),. {I = 7/2 S = 1/2) 

1_4 

when dissolved in liquid crystal media. These media inhibit the tumbling 

of the probe molecule to various extents depending upon the temperature and 
solvent material. If the probe is not completely free to tumble,’ then the epr 
spectra will be affected. The assumption that the spatial and temporal averages 
are equal breaks down, and the epr data must be interpreted in this light. 

The end cases of extremely fast and slow tumbling have been discussed in the 


literature 


3,4 


and are well understood. 


The intermediate case is more -difficult , 


however. Other authors^’ have successfully dealt with the intermediate case 
for isotropic liquids. The model which we are presenting is easily adapted 
to both nematic liquid crystal and isotropic solvents. The technique is also 
easier to deal, with from a computational point of view. 


In this model the molecules are assigned a discrete number of allowed 
angular positions. These positions are described by angles 6^ , 0^, etc. 

measured from the preferred direction in the case of a liquid crystal solvent 
and from the polar axis for the case of an isotropic liquid. Both axes cor- 
respond to the direction of applied magnetic field. A molecule in the 0^ 
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position can only relax to the 0^ ^ or the e n+ 2 _ position with a character- 
istic time x(0)o The choices of 8's and x’s are described below. The 
starting Hamiltonian in laboratory coordinates is given by 

# = ggH S + al °S + ^{AggH + bl ) ( 3 cos 2 8-l)S - ^<3 cos 2 6 -1 ) ( SI + s I ) 

( 1 ) 


Where g, a and Ag,b are the isotropic and anisotropic portions of the 
"g" and hyperfine tensors, respectively. S and. I refer to the electronic 
and nuclear spin operators. Equation (l) can be rewritten as 

#= aS + A(9) IS + f (SI + S I ) - B ( 0 ) (SI + S I ) (2) 

: z z z 2 + - -+ +- - + 


where 


a = 
B(0) 


[ 


1 + \ ( 3 cos 2 0-1) 

3 g 


Ijj (3 cos 2 0 - 1) 


gaH z’ 


and 

A( 0 ) = a + ^ b(3 cos 2 0-1) 

An rf field is applied 

H = u) (S cos ut - S sin wt) (3) 

_L -L x y 

Using the density matrix notation in the rotating frame we have 

p r = -nTh^te), p r ] (k) 

where 

*£-*&♦* (5 > 


and 

= (o-u)S z + A(0)S z I z - ul z + |<S + I_ + S_I + ) - B(0)(S + I_ + S_I + ),. 

H 1 = ML S x (^) 

Let p = p + p, , where p is zero order in to, and p, is linear in to, . 
r o 1 o 11 1 

In addition for a steady state solution p^ =0. Substituting into (L) and 

retaining terms up to first order in 
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p 1 (0) 

A term of the form — must be included to. represent relaxation of the 

2 1 

molecule. due to random perturbations and a term, - to represent relaxa- 

tion to other angles. The form of t( 0) is obtained using a diffusion approach, 
The 0 component of the diffusion equation for the isotropic solvent' is 


(dv 2 v) 


D 


0 component 


9 

30 


r sin 0 

where V is an arbitrary scaler field. Let 


— — 

. „ 3V 

D 

av a a v ~| 

sin 6 30 

* 2 
r 

L 8 « + 


( 8 ) 


A = 0 . - 0 . 

x+1 1 


then 


av v(e. +1 ) - vQ,-,) 

30 2A 


(9) 


) 

O 


and 


3 2 v T(6 i+1 ) -27(8^ * V(9 1 _ 1 ) 


a a a 

so that equation (8) becomes 


(dv 2 v) 


0 component 


1 

T 


2V(0 i ) --V(0 1+1 )< 1 + 2 cot 6 


where 


D 


t 2 2 
r A 


Using the format of (ll) 


- V( 0 


, )f i_L 

i-l ; \ 2 


*} 


cot 0 . 
i 


( 10 ) 


(11) 


xXeTT x 


2p(0 t ) - p(e i+1 ) ( i + g cot e i 


- p(. 0._ 1 ) i - f cot 0. 


( 12 ) 







Equation (7) is now modified to become, for the i — position, 
,-l 


0 = -iK 


^o (6 i )>p l (e i ) 1 ^d. (6 i ) » p o (e i ) 


- P,(9,) 


i (0 i ) 1 r 

T 2 " t 2 p l^ 9 i^ _ Pi 


(0 i+ l ) 1 + 2 COt 6 i 


- p l (e i-l } 1 -2 COt 6 i 


(13) 


0^,P q ] can be evaluated as 


. co. u 
1 1 o 


[S + - Sj, 


2 kT 

where co = g£5H . Equation (13) now becomes 

r I i P i < e i ) i f 

j^iWVJ -^--7 L 2p i (e i ) " p i (e i 


(14) 


i+1 


1 + — cot 0 . 

2 l 


~p 1 (e - i) 


1 - | cot 




(15) 


where sty is now in radians sec 

'o 


-1 


It is now necessary to choose values of 0 . In order to a.void computa- 
tional difficulties, 0 = 0, i.e., cot 0 = °°, will be avoided. We choose 
0^ = tt/ 34 and A = tt/1T. This gives 9 allowed positions in the first quadrant 

with 0^ = ir/2. Since the equations are independent of <f>, p(0^) = p(0^)„ We 

also assume that the molecules cannot be distinguished end from end and this 
yields p(0.) = p(ir + 0.). This makes the four quadrants equivalent so that 
we will restrict ourselves to the first quadrant. Equation (15) then is really 
a series of 9 equations corresponding to the 9 allowed positions. 


The absorption is proportional to 
M = Tr p 1 (0 1 ) S = - |- Tr Pl (S + - S_) . 


(16) 


VAAC has a nuclear spin of 7/2, thus its epr spectrum consists of eight lines. 
The line corresponding to = J/2 is obtained as 


absorption = Img <7/2, 1/2 | p (0 ) j 7/2 , -l/2>sin 6 ^ 


(17) 
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The forms for the other values of m^'s are similar. The m_j. = 7/2' s line 
is coupled to the m^. = 5/2's line via the I operator so that it is necessary , 
to. solve an 18 x 18 matrix to obtain the <7/2 ,1/2 1 p ( 0^ ) 1 7/2 ,. - l/2> .elements. 

A 36 x 36 matrix is required for all of the other lines except nu = -7/2 
since these lines couple to the higher and lower hyperfine lines via the I 

operators. T ? is obtained from the experimental line width,' and x is varied 

^ -12 -7 
from approximately ,10 sec for the free tumbling case to approximately 10 

sec for the slow tumbling case. An attempt will :be made to correlate the 

diffusion constant D with the viscosity via the Einstein relation • 


D = 


kT 


6iTari 

where n is the viscosity and a is the effective radius. 


(18) 


Three changes must be made to apply this model to a nematic liquid crystal <= 

First, in equation (17) the <m_,l/2jp(0. ) |m T ,-l/2> summation must be multiplied 

by the Saupe distribution function exp(-qcos 0/kT) as well as the polar 
volume element. Second, the form of the relaxation time must be changed to 
reflect the fact that the viscosities are different in the parallel and per- 
pendicular directions. Once again, taking the diffusion approach 



+ D„ 



and proceeding in a manner similar to the -isotropic case one obtains 


(19) 





1 _ 1 
t( 0) x 


D 2- D l . 2 


2p(eJ{ 1 + — - isin" 6, \ - p(e,^)U + 


V D 1 . 2 


i+1 


D., 


sin 0 . + — cot 0 „ 
1 2 x 


+ — 


A D 2" D 1 . 


2 D. 


sin 20^ > - p(0^ ( 1 + 


D 2- D l . 

D, 


sin 0. - — cot 0. 
1 2 1 


A D 2", D 1 . „ 

2 -X ~ Sln 29 i 


J 


( 20 ) 


This reduces to the isotropic case for Third, x(0) must be multi- 

plied by the Saupe function to, reflect the different populations at the different 
angular positions. 


The epr. spectra will be studied as a function of both the order parameter 



and. the relaxation time, and an effort will be made to separate the effects 
of these two parameters on the data. We will also. attempt to correlate the 
diffusion constants . and with the perpendicular and parallel components 

of the viscosity via the Einstein relation. 
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